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Abstract 

A popular model of decoherence based on the linear coupling to harmonic oscillator heat baths 
is analized and shown to be inappropriate in the regime where decoherence dominates over energy 
dissipation, called pure decoherence regime. The similar mechanism essentially related to the 
energy conservation implies that, on the contrary to the recent conjectures, chaotic environments 
can be less efficient decoherers than regular ones. Finally, the elastic scattering mechanism is 
advocated as the simplest source of pure decoherence. 
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Decoherence became one of the most popular topics in the physical literature of the last 
decade [1-4]. This is mainly due to the progress in experimental techniques allowing to 
observe the onset of decoherence at the most interesting regime i.e. at the border between 
quantum and classical worlds [5]. Another motivation is a destructive role of decoherence in 
the possible future technology based on quantum information processing [6]. Despite the fact 
that the theoretical models of decoherence exist at least for 40 years [7] a closer look at certain 
aspects of these theories reveals inconsistencies and interpretational problems. We shall 
use the following definition: Decoherence is an irreversible, uncontrollable and persistent 
formation of quantum correlations (entanglement) of the system with its environment. 

Usually, decoherence is accompanied by dissipation i.e. the net exchange of energy with 
environment. For the sake of clarity we shall restrict ourselves to the case of pure decoherence 
called also dephasing for which the process of energy dissipation is neglible. This situation 
occurs in two cases: A) for system's Hamiltonian Hs commuting (approximatively) with 
the system-bath interaction Hamiltonian Hint] B) for the initial states of S which evolve 
very slowly under the dynamics governed by Hg on the time scale of decoherence processes. 
In both cases we can disregard the presence of system's Hamiltonian Hs in the derivations 
concerning decoherence processes. As an example of B) we can consider a system equivalent 
to a 1-dimensional particle in a symmetric double-well potential coupled to a bath. In the 
semiclassical regime we can rectrict ourselves to the 2-dimensional Hilbert space of initial 
states spanned by the lowest lying almost degenerated Hamiltonian eigenstates [8]. The 
second example is a heavy particle interacting with a medium. As initial states for its center 
of motion we can choose superpositions of well-localized wave packets with small enough 
kinetic energies. For decoherence effects grow with the size of the particle while the kinetic 
energy exchange decreases we are again in the pure decoherence regime. 

Pure decoherence is supposed to be the main ingredient of the theory explaining the ap- 
parent absence of superpositions of macroscopically distinguishable states and the transition 
from quantum to classical world. Indeed, the explanation of the rapid decay of quantum 
correlations between "Schrodinger cat" states should not essentially depend on the energy 
difference between \dead cat > and \alive cat > but is rather related to the distinguishability 
of these states described in terms of certain collective observables which are coupled to the 
environment. 

In this letter we show that the most popular model of quantum open system based on 
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the linear coupling to the harmonic oscillator bath is inadequate in the pure decoherence 
regime. Using a unifying approach in terms of reservoir's spectral function we also show 
that, contrary to the recent conjecture, chaotic environments can be less efficient decoherers 
than regular ones at least for pure decoherence case. 

The first problem can be understood in simple physical terms. Namely, pure decoherence 
in the open system must be accompanied by the irreversible perturbation of the environ- 
ment's state but the energy of the environment should be asymptotically preserved. How- 
ever, the linear coupling to the bosonic environment implies that the only change of its state 
is caused by irreversible processes of emission and absorption of single bosons which must 
alter the environment's energy. The energy exchange can be reduced by a strong coupling to 
low energy bosons what in turn produces infrared divergencies. Those divergencies change 
completely the physical interpretation of the model, in particular the decomposition of the 
total system into the open system S and the reservoir R, what seems not to be taken into 
account in the literature [9,10]. 

Spin-boson model [11]. To explain this effect in a rigorous way we consider a two level 
system coupled linearly to the bosonic reservoir that is defined in terms of fields satisfying 
canonical commutation relations [a(a;), 0^(0;')] = S{u! — ui'), a single-boson Hilbert space 
L^[0,oo), a single-boson Hamiltonian hi, {hif){uj) — ujf{u!) and the second quantization 
Hamiltonian 

roo 

Hb — / duj uj aUu)a(uj) (1) 
Jo 

acting on the bosonic Fock space JFb(L^[0, 00)) with the vacuum state Q. The general spin- 
boson Hamiltonian depending on the function ("formf actor") g^cu) can be written as {a^ - 
Pauli matrices) 

Ug ^ Hs + Hb + Hintig) : Hs = \eai (2) 
where the interaction Hamiltonian is linear in bosonic field and reads 

Hint{g) = (73 / dwu{g{ijo)a{u) + g{u)a}{u)) . (3) 
Jo 

The Hamiltonians act on the Hilbert space TIsb = ^b(L^[0, 00)). The system Hamil- 
tonian Hs describes the coherent tunneling between two eigenstates of the " position" oper- 
ator (T3'0± = iV'i- We assume that this process is much slower then the decoherence i.e. 
h/e» Tdec ■ This is precisely the condition for pure decoherence regime in our model which 
allows to put e = in (3) and leads to 
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H, = ( ' ^ J ^ diag[i/+„ (4) 



where 



TOO I'OO 

H±g— I duj uj aUijj)aiijj) ± / duuigUjj)aUjj) -\- gUjj)aUuj)) (5) 

are van Hove Hamiltonians for the bosonic field which are well-known exactly solvable toy 
models of renormalization, both in the infrared and ultraviolet regimes. As pure decoherence 
is a low energy phenomenon we introduce an ultraviolet cut-off {giuj) = for a; > uJc) and 
for simplicity we assume a power-like low energy scaling 

\g{uj)\^ r-^ ijj'''^ ,ioY uj « UJc ■ (6) 

The most frequently used in the context of decoherence is the ohmic coupling i.e. k — 0. 
Under the assumptions of above one obtains the following rigorous results concerning the 
existence and properties of the van Hove Hamiltonians [12]: 

1) K > 0, regular case, ground states for H±g exist, 

2) —1<K<0, infrared problem, (includes ohmic case!), H±g - bounded from below but 
ground states do not exist (in the Fock space) , 

3) K < —1, uphysical case, H±g - unbounded from below or do not exist as self-adjoint 
operators on the Fock space. 

The main tool used in the analysis of the Hamiltonian is its diagonalization in terms of uni- 
tary Weyl operators W{f) — exp{a(/) — «'*'(/)} with ||/||^ = /o°° \ f{u!)\'^du < oo , acting on 
the Fock space and satisfying Weyl commutation relations: W{f)'^ — W{—f), W{f)W{h) — 
e-'^'^<f'''>W(f + h), W(f)a(u;)W(f)^ = a(u;) + f(u;)l. The vectors W(f)n are called coher- 
ent states and they form an overcomplete set in the following sense. If for a given vector ^ 
from the bosonic Fock space and any / e L^[0, oo) < \E', W{f)fl >= 0, then ^ = 0. Taking 
into account Weyl relations and < fl, W{f)fl >= exp{ — (l/2)||/p} we obtain 

I < W{f)n, W{h)Q > |2 = e-ll^-'^ll' < e-(ll^ll-ll'^ll)' , (7) 

so W{f) is not unitary on the Fock space unless ||/|| < oo. Introducing now for a given form- 
factor g{uj), \\g\\ < oo the unitary operator on Hsb defined as W{g) — diag[l^(g'), 1^(— g')] 
we obtain the diagonalized form 

W{g)UgW{gy ^ dwua\u)a{u)-Eg (8) 



where Eg —< g,hig >— ijj\g{Lo)\'^duj . Therefore, the degenerated ground states of 
are given by 

H,$±(^) = -^^$±(5) , ^±{g) = 'il^±®W{±g)n . (9) 

Obviously, in our setting the regular case > 0) corresponds to the condition \\g\\ < oo. 
Two degenerated ground states of the Hamiltonian Hg should be interpreted as the states 
of a dressed spin which consists of a bare spin and a cloud of virtual bosons represented by 
the coherent states W{±g)Q. The dressed system is decoupled from the environment and 
its dynamics is trivial. 

The case 2) (—1 < k < 0) corresponds to \\g\\ = oo, Eg < oo. In principle, the states 
^±{g) treated as limits of "normal" states from the Hilbert space Hsb with ||gf|| oo can 
exist in the sense of state functionals on the algebra of observables. But in this case they 
are disjoint, i.e. they define nonequivalent representations of the algebra of observables (van 
Hove phenomenon [12]). Formally, it follows from the formula ^+{g) = Ui (8) W{2g)^_{g) 
which for \\g\\ oo indicates that there exists no unitary operator which transforms ^-{g) 
into [13]. Physically, it means that their superpositions are indistinguishable from 

their mixtures {supers election rule). Some authors invoke this mechanism to describe the 
emergence of classical observables for quantum systems [14,16]. This phenomenon should 
be called static decoherence because the disjointness is a permanent feature of these states. 
Although from the mathematical point of view this is an atr active approach, on the other 
hand it can lead to profound interpretational difficulties, e.g. unphysical superselection rules 
in quantum electrodynamics[16,17]. 

We invoke now the standard dynamical approach to decoherence in open systems applied 
to our spin-boson model. As an initial state of the total system one chooses the product 
state 

= ^ (g) Q , ^ = a^ijj_ + q;+^+ , a± eC (10) 

satisfying 

I < *m|$±(5) > I' = e-ll^ll' , E{^in) =< *m|H,|*,„ >= , (11) 

computes its time evolution governed by the Hamiltonian (4) (5) 

^(f) = e-'^^^^in = exp{i{tEg - Im < g\gt >)}(a_V'- W{gt - g)^ + a+i^+ (g) W{g - gt)^) 

(12) 
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where gt{uj) — e '^'^*g{uj) and calculates the reduced density matrix for spin 

p, = IV5|*(i)><*(i)|= _ (13) 

\a+a_e '* lo^-l / 

with 

7t = 2||^-^t|p<8||^|p. (14) 

Usually, one discusses the structure of the reduced density matrix (13) only, and identifies 
74 with the decoherence factor. It follows from (14) that decoherence is complete only if 
llg^ll = oo i.e. for a singular coupling. To obtain an asymptotically exponential decay of the 
off-diagonal elements of the reduced density matrix (13) we must assume 

0<7=lim^=lim ru;2| (u;)|2i^^^da; = 7rlima;2| (^^|2_ ^^^^ 

This result agrees with the standard wisdom relating the pure decoherence rate to the value 
at a; = of the spectral density function [4,18] 

R{u) - / e*'"* < R{t)R >B dt (16) 

where i? is a bath's operator appearing in the interaction Hamiltonian a^® R and < ■ >b 
is an average with respect to the environment's state. It is a special case of the quantum 
fluctuation - dissipation theorem which in fact should be called in this context a " fluctuation- 
decoherence theorem". For our model R — J du ijj{g{uj)a{uj) + h.c.) and hence Ro{u) — 
27ra;^|gr(a;)P where the subscript "0" indicates the zero-temperature (vacuum) state of the 
bath. However, a non zero value of 7 means k — —\ which is the uphysical (subohmic) case 
2). The situation is slightly less singular for the temperature T > 0. It is easy to check 
that for u « T (we put h = ks = 1) -Rt('^) — {T/u!)Ro{uj) and hence the condition of 
finite decoherence rate is satisfied for the ohmic case k — 0. The same ohmic assumption is 
made in the derivation of the popular Caldeira-Leggett equation for the quantum Brownian 
particle coupled to harmonic oscillator heat bath [10]. 

A different physical interpretation of the discussed process follows from the analysis of the 
state evolution for the total system in the regular case. As for t 00 the traveling wave 
becomes orthogonal to g the asymptotic form of ^(t) possesses the structure of superposition 
of two triple product states ip± W{±g)fl W{^gt)fl. Hence, the evolution of the initial 
product state (10) given by (12) describes the process of formation of the cloud accompanied 
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by emission of the average energy Eg in a form of coherent travehng waves ±gt. Therefore, 
from the physical point of view the discussed model describes a phenomenon which should 
be called false decoherence [19]. 

Summarizing the above results one can say that the linear coupling to harmonic oscillator 
bath is not an appropriate model of pure decoherence phenomena. In general, we have two 
different choices of the formf actor g{u!) - a singular and a regular one, both leading to 
difficulties. 

The singular couphng (e.g. ohmic, k = 0, for T > or subohmic, k, — —1, for T — 0), 
which gives a finite asymptotic decoherence rate in the formal derivation of the reduced 
density matrix, leads to serious interpretational difficulties due to infrared divergencies. The 
ground states do not exist in the Fock space, those defined by limiting procedures describe 
disjoint states for which superposition principle is not applicable. As a consequence the 
picture of an open system gradually loosing its "quantum coherences" is not valid in this 
case. 

For the regular coupling the natural representation of dressed states gives the picture of a 
"physical dressed system" decoupled from its environment. The choice of the initial product 
state of a bare system and a bath makes sense only in the unique moment of system's creation 
followed by the irreversible dressing process. After that, the system cannot be prepared in 
a such state again. The absence of true decoherence in the regular (superohmic) case is 
indicated by the relation limi^^o -R(<^) = 0. It follows that the models based on vacuum 
fluctuations of the background quantum fields (gravitational, electromagnetic,...) [20] are 
very unlikely to solve the problem of transition from the quantum to classical world. 

Chaotic vs. regular environments. The similar physical mechanism leading to the ab- 
sence of pure decoherence appears in the case of a bath which is an ensemble of quantum 
subsystems with chaotic properties. The intuition supported by some heuristic arguments 
suggests, as it is formulated in [21], that "...one would expect that environments with un- 
stable dynamics will be much more efficient decoherers,...". However, it follows from the 
previous discussion that the reservoir's energy eigenstates should be degenerated and la- 
beled by other quantum numbers which can be altered without energy modification. On 
the contrary, for a chaotic system its energy levels are typically nondegenerated due to the 
mechanism of level repulsion [22]. 

We begin with the physical example as a motivation. Consider a large system (say a 
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molecule) with a relevant collective, single degree of freedom which can be modeled by a 
2-level system as above. This degree of freedom is our open system again, while the internal 
degrees of freedom of the molecule form a bath which can be treated as a large ensemble of 

quantum systems. 

The simplified mathematical model consists of a 1/2-spin system interacting with an 
ensemble of identical M-level quantum systems by means of the following mean-field 
type Hamiltonian which is an anolog of (2) [23] 

1 N N 

Hq = ^ea, + Y: h^^^ + <73 ® N-'l^ Q'^'^ ■ (17) 

^ k=l k=l 

Here /i*^*^^ is a copy of the Hamiltonian with the spectral resolution h — J2m=i ^ml"^ >< "^|, 
> em and Q^^) is a copy of an operator Q — Q\ TrQ — 0. The rerefence state 
of environment is assumed to be the product state where p^'^^ is a copy of the 

microcanonical state giving an uniform probability distribution over all states \m >. We 
assume again, that the tunneling time h/e is much longer that the dccohcrence time and 
that the energy e is much smaller than the typical reservoir's constituents energy spacing 
A (compare again [8]). For N ^ oo the mean-field reservoir's observable N~^I'^Y,^^-^Q^^^ 
behaves like a Gaussian noise and in the Markovian approximation the pure decoherence 
rate 7 for the spin is given by the following version of the fiuctuation-dissipation formula 

1 - - TT ^ 

J = -lim R{uj) , R{u}) ^ — Yl \<m\Q\m' >\'^5{{em-em')-ou) . (18) 

m,m'=l 

This formula makes sense also when instead of identical subsystems the reservoir consists of 
a large random ensemble of quantum systems with Hamiltonians h^''^ characterized by the 
average nearest-neighbour level spacing distribution p{s) with the average A. For u << A 
only the nearest-neighbour level spacings s = tm+i — Cm contribute to the spectral function 
R{u!) (18)and therefore the difference between the bath consisting of classically integrable or 
chaotic systems becomes crucial. For the former we expect a Poisson distribution of s while 
for the later the level repulsion given by p{s) ~ s'^, /3 = 1, 2, or 4 is generically observed [22]. 
Assuming that the magnitude of the matrix elements | < m + l\Q\m > \ is not strongly 
correlated with e^+i — we obtain 

R(u;) ~ 7rQ^p{uj) (19) 
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where is an averaged value of| <m + l|(5|m> p. As a consequence pure decoherence 
rate is equal to zero for the chaotic systems while for the regular ones we obtain a finite 
value of 7. 

We have shown that for important models used to analyse environmental decoherence 
this phenomenon disappears in the limit of pure decoherence. These cases can be easily 
detected applying the unified approach of spectral density function (16) at a; ^ O.To avoid 
the possible influence of approximation procedures we have studied in details an exactly 
solvable spin-boson model which illustrated both , vanishing pure decoherence and infrared 
divergencies for the singular choice of the coupling. The proper models of pure decoherence 
involve elastic scattering processes which "for all practical purposes " provide strong enough 
suppression of quantum superpositions in the macroworld (e.g. Joos and Zeh [7]), [24] or 
describe quantum Brownian motion [25]. 
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